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An adaptive spline is used in this work to deal with singularly perturbed 
boundary value problems with layers in the interior region. To evaluate 
the layer behavior in the solution, a different technique on a uniform mesh 
is designed by replacing the first-order derivatives with nonstandard differ- 
ences in the adaptive cubic spline. A tridiagonal solver is used to solve the 
tridiagonal system of the difference scheme. The fourth-order convergence 
of the approach is established. The validity of the suggested computa- 
tional method is demonstrated through numerical experiments, which are 
compared to other methods in the literature. Layer profile is depicted in 
graphs. 


AMS subject classifications (2020): 65L10; 65L1; 65L12. 


Keywords: Singular perturbation problem; Interior layer; Adaptive spline; 
Tridiagonal system. 


* Corresponding author 
Received 23 November 2021; revised 29 December 2021; accepted 29 December 2021 


Erla Srinivas 


Depar 


ment 


Hyderabad. 


Mudavath Lalu 


Depar 


men 


Hyderabad. 


Kolloju Phaneendra 


Depar 


men 


of Ma: 
email: 


of Ma: 
email: 


of Ma: 


thema 


hema: 


hema: 


ics, University College of Engineering, Osmania U 
srinivas.erlaQ@gmail.com 


ics, University College of Engineering, Osmania U 
lalunaikmud 


athou@gmail.com 


ics, University College of Engineering, Osmania U 


Hyderabad. email: kollojuphaneendra@yahoo.co.in 


355 


niversity, 


niversity, 


niversity, 


356 Srinivas, Lalu and Phaneendra 


1 Introduction 


A small parameter leads to singular perturbation problems (SPPs) in a va- 
riety of science and engineering problems in fluid mechanics, elasticity, aero- 
dynamics, magneto-hydrodynamics, optimal control, and other domains of 
fluid motion [1]. WKB difficulties, the modeling of steady and unsteady vis- 
cous flow problems having big Reynolds numbers, magneto-hydrodynamics 
duct problems with high Hartman numbers, and so on are few notable ex- 
amples. The numerical solutions to SPPs are often nontrivial because of the 
rapid development or decay (boundary/interior layer behavior) of solutions 
to dissipative problems. Here, we include a few references related to the in- 
terior layer problems. The authors in [2, 3, 5, 11, 16, 24, 25] provided a full 
theoretical, analytical, and numerical discussion of this topic. The authors 
of [8, 9] provided an informative overview of SPPs on layers in their survey 
publications. In [10], for the Emden—Fowler type equations, the author pro- 
posed a finite element collocation technique based on cubic B-splines. Farrell 
[4] provides a uniformly convergent scheme for the SPP with a turning point 
problem. Geng, Qian, and Li [6] suggested an approach based on the re- 
producing kernel method and the asymptotic expansion technique. Miller, 
O’Riordan, and Shishkin [12] explained how to solve convection-diffusion 
and reaction-diffusion problems using conventional techniques on Shishkin 
meshes. Natesan, Vigo-Aguiar and Ramanujam [13] split the region into two 
subdomains, the layer domain and the regular domain, and then resolved the 
layer problem with the fitted operator method and the regular problem with 
the finite difference technique. Navnit [14] proposed a fourth-order numerical 
approach for solving an SPP using an adaptive cubic spline. The same au- 
thor in [15] proposed a general scheme for the solution of nonlinear SPP using 
an adaptive cubic spline. The author in [19] developed a nonuniform mesh 
optimal B-spline collocation method for the numerical solution of a singular 
two-point boundary value problem describing electro hydrodynamic flow of a 
fluid in a circular cylindrical conduit. In [20], two B-spline collocation meth- 
ods were proposed to solve a class of nonlinear derivative dependent singular 
boundary value problems. Ramos [22] presented a locally analytical tech- 
nique for solving SPPs with internal and boundary layers, as well as turning 
points. 

With this motivation, in this paper, a higher order finite difference method 
on a uniform mesh for solving SPP turning point problems with an interior 
layer is proposed. 


2 Statement of the problem 


Consider a second-order SPP of the form 


ew’ (t) + P(t)w'(t) + Q(t)w(t) = fd), -1<t<1l, (1) 
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with boundary conditions 
w(-l) =a and w(1) = 8, (2) 


where 0 < € < 1 is a positive perturbation parameter, a and £ are finite 
constants, and P(t), Q(t), and f(t) are considered to be suitably smooth 
functions. Based on the coefficient P(t), the solution to (2) has a layer or 
turning point behavior. The presence of a turning point in the solution to 
the problem makes it substantially more difficult to handle a boundary or 
inner layer. In this paper, we examine at the situation where a problem’s 
turning point results in a solution with an interior layer. Under assumptions 


P(0) =0,P'(0) 20, Q(t) <Q <9, for allt € D = ([-1, 1, 


for allt € D = [-1, 1], 


the given turning point problem has a solution with interior layers at t = 0. 


3 Adaptive spline 


With grid points t; in [a,b], consider a mesh such that {Q: a = tp < th < 
to <-++<ty =b, where h = t;—t;-1 fori =1,2,...,N}. A function w(t,7) 
interpolates w(t) at the grid points t;, depends on a variable, and leads to a 
cubic spline w(t) in [a, b] as 7 + 0, named as an adaptive spline; see [7, 23]. 
The function w(t, 7) satisfies the equation 


el" (t,7) -pW(i7) = 


t; —t 
(eM; — pmi)+ je (6¢Mi-4 = pig) (3) 
where f;-1 <t < tj, w’(t,7) = mj, and w’(t,7) = M;. Solving (3) and using 


the interpolatory conditions w(t;-1,7) = w;-1 and W(t;, T) = wi, we have 
w(s,7) = A; + Bye™* — Ke [s7?22 +7241] (Mi — F<mi) 
h? [1 2 3 T 
Ta 57 (lL=@ sels) el (Min - rmi-1) ; (4) 


where 


A; (e7 — 1) = —x; + aj_-1e —3 (S+7+1) —Te | (Mi - Emi) 
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h2 
B;(e7 — 1) =a; — aj_1e7 + 3 lg +1) —Te™| (M; - - mi) 


2 h 
+ G aa 1) = 7] (Mi-1 = ;mi-1), 
— Phe and z = *—'+. The spline function ~(¢,7) on [t;, ti41] is acquired 


with replacing i by (¢+ 1) in (11). Utilizing the first or second derivative 
continuity conditions of a(t, 7) at t = t;, we get the following relationship: 


T (i 
(Mins = omit) los +7r4+1)- q 


es [e~ (5-7-2) + (-5-++2)| 


+ (Min = imi) le" —l+rTr-— 5 


2 
= aie [e-"wiga — (1+e77) wi; + wi-a]- (5) 


Furthermore, the below relations are given for the adaptive splines: 


1 
(7) mi—-1 = —h (Ar Mj-1 + A2Mi) + i (w; — Wi-1), 
1 
(3) my = h(A3Mi-1 + AaMi) + i (wi — wi-1), 
oh i — Wi- 
(iii) M4 Sana Pain) p ain 1) 
2T h 
Oh 4 Wie 
(iv) 5 Mi = (Asmi-1 + Armi) + Bo Mee W) 
2T h 
where ; ; : , 
Aca Girt Aaa os 
1 6 1 6 
Ag=7(l+8)-s, A= 7-9) +5; 
B =+(1-6) B oe 2G 429) and 9 =coth(5) — = 
1 4 ; Dire 2 ’ 1 = co 5 = 


We also obtain 


1 
A2Mi41 + (Ai + A4) M; + A3sMj-1 = h2 [wi4d 2w; 4 wi-1| . (6) 
Remark: In the limiting case when 7 — 0, we have 


1 1 1 1 0 
Ai = As= 3; A2z=A3=5, Bi=5 = 


and (8) reduces to the ordinary cubic spline scheme. 
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4 Description of the procedure 


At the mesh point t;, the suggested approach can be discretized by the 
convection-diffusion equation (2) as 


eM; = f(ti) — P(ti)wi' (t) — Q(ti) w(t). (7) 
The above equations shall be replaced by (8), and by approximating the first 
order derivatives of t at the mesh points t),to,...,ty_— 1 as 
1 —wi41 + 4u; — 8uj;-1 / 3wi41 — 4wj t+ wWi-1 
Wi-1 © oh ) Wi41 oh , 
1+ 2h? Qiia + 7h [3Pig1 + Pi 
wh ( nh? Qi+1 oA [BPis4 MD) wigs — 20 [Pig + Pia] wi 
1+ 2nh?Q;_1 — nhl P41 + 3B 
( a ie m7 Pit 1) wiit7 h[fisa — fia, 


we get the tridiagonal system 
Dywj-1 + Cj wi + Ujwi41 = Wi for i=1,2,..., N-1, (8) 
where 
L,=- e~ 5g Pith 
— (Ar + Ag) Pik [1 + 29h? Qi-1 — nh (Pig + 3Pi-1)] 


A 
+ ay tie + As Qi-1h? 


C; =2e + 2A3 Pi-1 h-4 (Ai + A4) Pih?n[Pi41 + Pi_4| 
— 2A, Pi4ih +2 (Ai + A,) Q;h? 


A 
U; =—-eE—- a P;_1 h + (Ai + A,) P;h [1 + 2nh?Qi-1 + nh (3Pi41 + F249) 


3 
+ gAanFitih + ArQizih? 


W, = h? [(A2 — 29 (Ai + Aa) Pih) fig +2(Ar+ Aa) Si] 
+h? (As + 27 (Ai + Aa) Ph) fi-i]. 


The tridiagonal system (10) is solved for i = 1,2,...,N —1, to obtain the 
approximations w1,w2,...,Wn_—1 of the solution w(t) at ti, te,...,tnw—1. 
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5 Truncation error 


Developed local truncation error associated with the scheme in (10) is 


T; (h) =e [1 — (2(Ar + As) + An + As)] A ad + ¢€ (As — Ao) how!” 


+|2 ase — Ane (Ai + As) -— 5 = (A + Aa) + Az + As]| Pirw”” 
Be [1 — 6 (Ap + As)] how?” 
1 

— 5 (As — Aa) [Pie +2(P! + Qi) wl” + 6(P! + Q)) w 


42 ne ni 3Q)”) ww +4 [2Q/”" we asa h? 4 O (n°) . 


Thus for different values of Aj, A3,A1 + Ag in the scheme, (10) indicates 
different orders. 

Remarks: 

(i) If Az = As, for any choice of arbitrary A2,.A,+A4 with (Ar + Ay)+A2 = 
$ and for any value of 7 , then the mewhod is obtained for second order. 

(ii) For Ap = A3 = 4, (A1+As) = 3, and n=- the fourth-order 
method is derived. 


one , 


6 Convergence 


The convergence analysis of the proposed method is discussed in this section. 
The system of equations (10) in the matrix form with the boundary conditions 
is 


(D+ F)W+G+T(h) =0, (9) 
where 
2e -- 0 0... O 
—e2-e 0... O 
Gee. Ce cee 0 -—e 26 -e... 0 
O> dienes 0 —e 2e 
and 
v, Ww, 0 0... 0 
Z Vo We 0 0 
Pana ats 0 23 U3 W3 ... 0 
Ore sacteee 0 ZN—1 UN-1 
in which 
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ss 3 
A= 74s P43 h = (Ai + A4) Ph [1 + Qnh?Qi-1 a nh (Pi41 + 3P;-1)| 


A 
+ oa 41h + Az Qi-1h? 


Uj =2A3 P;_4 h—-A4 (Ai + A4) P;h?n Lei +. P7| = 2Az Pyih 
+2 (Ai + A4) Q;h? 


s A 
Wi=— = Py h + (Ai + Aa) Pi [1 + 2nh?Qi-1 + nh (3Pi41 + Pi-1)] Pi 


a SAaPisth + AsQisih? forall 1=1,2,3,4,...,N—-1, 


and 


G= [1 = Z1Q, 42,93, +++)QN-1 — Wn-183] ’ 
in which 


Gi =h? [(Az — 2n (Ar + Aa) Pih) figa +2(Ar + As) fil 


+h? [(A3 + 2n (Ai + Aa) Ph) fi-1] for alli = 2,...,N—1, 
a: ae eee egg tia seh 
T(h) =0(h"), Ag = A3 = 12’ (Ay + Ag) = ig = He 


W=(Wi. Wa, Ws,-..,Weeal, Th) = (1, Ty, Trea , 
O =(0,0,...,0]". 
Let w = [wi, We,... weal” ~W satisfy the equation 
(D+F)w+G=0. (10) 


Let E = [e1,€2,.. w@neail = w—W, where e; = w;-W;, 1 = 1,2,...,N—-1, 
is the discretization error. Using (9) and (10), we obtain 


(D+F)E=T(h). (11) 


Let |Q(s)| < € and |P(s)| < &, where €; and & are positive constants. If 
¢i,j is the (i, 7)th element of F’, then 


[Gita] = |Wil < © (h (Ag + (Ar + Aa) €1 + hr? An€s + 4 (Ai + Aa) 0h? €7) 
+ (2h? (Ay + As) n€1€2) for alli =1,2,...,N —2, 

IGie—1] = [Zi] < © (h(Az + (Ai + Aa))m +h? Anne + 4 (Ai + Aa) nh?) 
+ (2h? (Ay + Ag) n€1€2) for all i = 2,3,...,N—1. 


Hence, for sufficiently small h, we have 


Cita <eé for all i = 1,2,...,N — 2, (12) 
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\Cii—a| < E for all i = 2,3,...,N— I. (13) 
Therefore, (D + F) is irreducible (see [27]). 
Let the sum of the values of the ith row of (D+ F) be S;. Then 
Aoh 2 
Si = oar a (Pyaa — 3P1)+h (Ai + Aa) P+h* (AoQi—-1 + 2 (Ar + Aa) Q:) 


+h? (Ay + Aa) nP; (3P)-1 + Pii1)— 2h? (Ar + Aa) nPiQi-1 fori=1, 
Si = h? (Qi-1 + 2 (Ar + Aa) Qi + AoQi41)+2h3 (Ai + As) Pin (Qi41 — Qi-1) 
for all i = 2,3,...,N—2, 
Ash 


Si =et+ 5 (Pi-1 — 3Pi41)—h (At + Ag) Py +h? (AoQi-1 + 2 (Ai + Aa) Qi) 
—h? (Ai + A4) nP; (3Pi41 + P,_-1) _ 2n3 (Ai + A.) nPiQi-1 
fori= N-—1. 


Let i+ = min, |P(ta)l let 7 = ymax, [Pal let 9+ = nin 1Q(ti)l and 
let 3 = max |Q(t)|- 

Since € is very small and € « o(h), for suitable small h, (D + F’) is mono- 
tone (see [26, 27]). Hence, (D+ F)~* exists and (D+ F)~* > 0. Thus from 


(14), we get 
ZI] < (D+F) NITIL- (14) 


Let (i, k)th element of (D+ F)~" be (D + F); ;, and define 


N-1 
—1y) —1 _ 
KD FY = gaa YD +P ice and (ION = gone POL 


Since (D+F); >Oand DY" (D+F); 4 5, =1, for alli=1,2,...,N— 
1, hence 


(D+ Pye < 5 < h? [(Az + 2(Ai + ie — 4(A, + Ag) v&?]’ eo) 
ew aoe s i] 
Sy —1 ~ hh? [(Ag + 2(Ai + Aa) bo — 4(Ar + Aad v2] 
Furthermore, for all 7 = 2,3,...,N — 2, . 
N-2 4 il 1 
k=2 ene ache 2 Sr < Be [2 (Az + (Ai + Aa)) &2*] os 


By utilizing (15)—(18) and using (14), we have 
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4 
|ZI| < O(*). (19) 


Hence, the method given in (10) is fourth-order convergent for 


1 5 
Az = As = 75) (Ai + Aa) = 55, and y=- — 


7 Numerical examples 


Three examples with internal layer behavior are examined in this part to 
explain the concept computationally. These problems were chosen because 
they have received considerable attention in the literature. The numerical 
rate of convergence is computed using the formula given by Doolan, Miller, 


(BE )_1lo 2N 
and Schilders [3] as. RN = 82" )les(2™"). 
Example 1. [17] Consider 
ew" (t) + 2tw’ (t) = 0, te (-1, 1), 


with w(—1) = —-1 and w(1) = 1. 
Thus w(t) = erf (+z) is the exact solution to the problem. 


Table 1 shows the computed solution’s maximum absolute errors (MAEs). 
Figure 1 graphically depicts the numerical and exact solutions. The error 
plot in the solution of this example is shown in Figure 4. 


Example 2. [21] Consider 
ew" (t) + 2(2t — 1)w’(t) — 4w(t) = 0, t € (0,1), 
with w(0) = 1 and w(1) = 1. Hence 


(1—2t)? (1—2t)? 


(224)? ) 
ewe Qe (2¢ ae V 2rter f (+) —e 2% w~W2nerf (42) a 2vé) 
ete /Inerf (Fe) + 2/e 


w(t) = 


is the exact solution to the problem. 

The MAE of the solution is shown in Tables 2 and 3. Figure 2 shows 
the layer profile of the numerical and exact solutions. The error plot in the 
solution of this example is shown in Figure 5. 


Example 3. [17] Consider 
ew’ + tw’ —w=0, -1<t<l, 
with w(—1) = 1 and w(1) = 2. 


The exact solution is 
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2 


2/e (t+ 3e~ "=" ) + (c% VOnterf (sz) 4 3e% V2ner f (+)) 
Qez= \/2rerf (Ar) +4 JE , 


w(t) = 


The MAE is shown in Table 4. Figure 3 shows the layer profile of the nu- 
merical and exact solutions. The error plot in the solution of this example is 
shown in Figure 6. 


Table 1: MAE of the solution of Example 1 


el N=2 N=25 N=2 N=28 N=29 N =210 
present method 
2-5 1.1881(-4 7.3536 (-6 4.5818(-7 2.8613(-8) 1.7889(-9 1.1198(-10) 
2-8 4.8395 (-4 2.9414(-5 1.8231(-6 1.1455(-7) 7.1558(-9 4.4732(-10) 
2-7 0.0019(-0 1.1881 (-4 7.3536(-6 4.5818(-7) 2.8613(-8 1.7890(-9 
2-8 0.0065(-0 4,8395(-4 2.9414(-5 1.8231 (-6) 1.1455(-7 7.1558(-9 
2-9 0.0203(-0 0.0019(-0 1.1881 (-4 7.3536 (-6) 4.5818(-7 2.8613(-8 
2-19 0.0957(-0) 0.0065(-0 4.8395(-4 2.9414(-5) 1.8231(-6 1.1455(-7 

RN 3.8757 3.7514 4.0403 4.0120 3.9923 

Results in [17] 
2-5 3.3962(-1 4.9291 (-1 6.1199(-1 7.0727(-1) 7.8166(-1 8.3818(-1 
2-8 2.6635(-1 4,2332(-1 5.5588(-1 6.6268(-1) 7.AT24(-1 8.1197(-1 
2-7 1.7078(-1 3.3962(-1 4.9291 (-1 6.1199(-1) 7.0727 (-1 7.8166(-1 
2-8 9.5121(-2 2.6635(-1 4,2332(-1 5.5588(-1) 6.6268(-1 7.AT2A(-1 
2-9 2.9280(-2 1.7078(-1 3.3962(-1 4.9291 (-1) 6.1199(-1 7.0727 (-1 
2-10 1.1065(-1) 9.5121 (-2 2.6635(-1 4,2332(-1) 5.5588(-1 6.6268(-1 


Total elapsed time with maximum number of subintervals 21° is 0.086497 sec. 


Table 2: MAE of the solution of Example 2 
el N=2 N = 26 N=2 N=2 N=2 N = 210 


present method 
2-5 8.6799(-6 5.4123(-7 3.3837(-8) 2.1151(-9) 1.3232(-10) 8.7138(- 
2-6 2.4578(-5 1.5331(-6 9.5748(-8) 5.9831(-9) 3.7396(-10) 2.3510(- 
2-7 6.7327(-5 4.3399(-6 2.7061(-7) 1.6919(-8) 1.0576(-9) 6.6184(- 
2-8 1.9559(-4 1.2289(-5 7.6654(-7) 4.7874(-8) 2.9916(-9) 1.8699(- 
2-9 5.6901(-4 3.3663(-5 2.1700(-6) 1.3531(-7) 8.4593(-9) 5.2878 (-1 
2-19 0.0015(-0) 9.7794(-5 6.1444(-6) 3.8327(-7) 2.3937(-8) 1.4958(-9 

R-N 3.9490 3.9924 4.0028 4.0010 4.0003 

Results in [18] 
2-5 5.9701(-3 3.3654(-3 1.7391(-3) 8.7449(-4) 4.3697 (-4) 2.1822(-4 
2-6 5.3525(-3 3.2322(-3 1.7219(-3) 8.7336 (-4) 4.3719(-4) 2.1834(-4 
2-7 1.1177(-2 2.9851(-3 1.6827(-3) 8.6953(-4) 4.3725 (-4) 2.1848 (-4 
2-8  2.5867(-2 2.6763(-3 1.6161(-3) 8.6093(-4) 4.3668 (-4) 2.1860(-4 
2-9  4.7842(-2 5.5886(-3 1.4925(-3) 8.4134(-4) 4.3477(-4) 2.1862(-4 
2-10 7.5829(-2) 1.2934(-2 1.3381(-3) 8.0805(-4) 4.3046 (-4) 2.1834(-4 


Total elapsed time with maximum number of subintervals 21° is 0.089720 sec. 
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Table 3: MAE of the solution of Example 2 
el N=% N=2 N=2" N=2 N=2? N = 210 
Present method 


2-© 2.4578(-5) 1.5331(-6) —9.5748(-8) —-5.9831(-9) —-3.7396(-10) —-2.3510(-11) 
2-8 1.9559(-4) 1.2289(-5) —-7.6654(-7) —4.7874(-8) 2.9916(-9) 1.8699(-10) 
2-10 0.0015(-3) —-9.7794(-5) —-6.1444(-6) —3.8327(-7) 2.3937 (-8) 1.4958(-9 
2-12) 4.4958(-3) ——-7.5528(-4) —- 4.8897(-5) —-3.0722(-6) 1.9163(-7) 1.1969(-8 
2-14 3.4263(-3) 1.7659(-3) 3.7764(-4) —-2.4449(-5) 1.5361/(-6) 9.5817(-8 
RN 0.9562 2.2253 3.9492 3.9924 4.0028 
Results in [21] 
2-6 3.630(- 3) 1.475(-3) 5.047(- 4) 1.508(- 4) 4.146(- 4) 1.089(- 4 
2-8 3.987/(- 3) 1.952(- 3) 9.168(- 4) 3.733(- 4) 1.272(- 4) 3.789(- 5 
2-19 7.147(- 3) 1.979(- 3) 9.814(- 4) 4.866(- 4) 2.297(- 4) 9.359(- 4 
2-  5.562(- 3) 3.583(- 3) 9.837(- 4) 4.898(- 4) 2.444(- 4) 1.215(- 4 
2-4 4,045(- 3) 2.778(- 3) 1.794(- 3) 4.908(- 4) 2.446(- 4) 1.222(- 4 


Table 4: MAE of the solution of Example 3 
el N=2 N=2 N=2 N=2 N=2 N= 210 
present method 


2-5 1,3020(-5) —-8.1184(-7) 5.0756(-8) 3.1727(-9) —1.9854(-10) —1.3311(-11) 
2-° 3.6866 (-5 2.2996 (-6 1.4362(-7) 8.9747(-9) 5.6092(-10) 3.5138(-11) 
2-7 1.0099(-4 6.5099(-6 4.0592(-7) 2.5378(-8) 1.5863(-09) 9.9321(-11) 
2-8 9.9338(-4 1.8433(-5) 1.1498(-6) 7.18118)  4.4873(-09) —_-2.8046(-10) 
2-9 8.5352(-4 5.0495(-5 3.2550(-6) 2.0296(-7) 1.2689(-08) 7.9318(-10) 
2-10 2.9657(-3)  1.4669(-4) —-9.2166(-6) 5.7490(-7) 3.5906(-08) —_2.2437(-09) 

RN 3.9491 3.9924 4.0029 4.0010 4.0003 

Results in [18] 

~5  1,0203(-2 5.5017(-3 2.7453(-3) 1.3449(-3) 6.6344(-4) 3.2928(-4 


2 2 
2-6 9.1104(-3 5.3503(-3 2.7665(-3) 1.3581(-3) 6.6733(-4) 3.3030(-4 
2-7 7.9110(-3 5.1014(-3 2.7508(-3) 1.3726(-3) 6.7246 (-4) 3.3172(-4 
2-8  2.2330(-2 4.5552(-3 2.6751(-3) 1.3832(-3) 6.7905(-4) 3.3366 (-4 
2-9 4.6794(-2 3.9555 (-3 2.5507(-3) 1.3754(-3) 6.8632(-4) 3.3623 (-4 
2-10 7.7601(-2) 1.1165(-2 2.2776(-3) 1.3376(-3) 6.9162(-4) 3.3953 (-4 
Results in [17] 
2-5 1.0111(-1 1.3778(-1 1.6169(-1) 1.7746(-1) 1.8802(-1) 1.9520(-1 
2-6 5,3820(-2) 8.5843(-2) —«1.0677(-1) —-1.2048(-1) —«1.2958(-1) ——1.3573(-1 
2-7 2.2863(-2 5.0556(-2 6.8889(-2) 8.0847(-2) 8.8730(-2) 9.4011(-2 
8 


5.6850(-3 2.6910(-2 4.2921(-2) 5.3386 (-2) 6.0238 (-2) 6.4792(-2 
2.0288 (-2 1.1431(-2 2.5278(-2) 3.4445 (-2) 4.0424 (-2) 4.4365 (-2 
2-10 2.9876 (-2) 2.8425(-3 1.3455(-2) 2.1461(-2) 2.6693(-2) 3.0119(-2 


Total elapsed time with maximum number of subintervals 2!° is 0.090109 
sec. 
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Figure 1: Numerical and exact solution of Example 1 for ¢ = 2719, h = 27-7. 
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Figure 2: Numerical and exact solution of Example 2 for e = 2719, h = 277. 
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Figure 3: Numerical and exact solution of Example 3 for ¢ = 271°, h = 277. 
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Figure 4: Log-Log scale for Example 1 
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Figure 5: Log-Log scale for Example 2 
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Figure 6: Log-Log scale for Example 3 
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8 Conclusion 


In this paper, for solving singularly perturbed two-point boundary value prob- 
lems with an interior layer, a higher-order finite difference approach was 
proposed. To derive the discretization equation, an adaptive cubic spline ap- 
proach has been extended for a singularly perturbed boundary value problem 
with interior layers. It was produced by substituting nonstandard finite dif- 
ferences in the first-order derivatives of w(t). The tridiagonal solver discrete 
consistent embedding was used to solve the discretization equation. Conver- 
gence was evaluated in the proposed method. To demonstrate the method, 
numerical problems have been solved, when ¢ is either small or large as com- 
pared to the mesh size h. To justify the method, numerical results were 
compared with the results taken by the methods given in [17, 18, 21]. For 
the layer behavior, the solutions were graphically displayed, and we discov- 
ered that the numerical solution closely matches the exact solution. This 
technique is simple to compute and requires little computational effort. 
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